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ABSTRACT 

The opening of DNA double strands is extremely 
relevant to several biological functions, such as rep- 
lication and transcription or binding of specific 
proteins. Such opening phenomenon is particularly 
sensitive to the aqueous solvent conditions in 
which the DNA molecule is dispersed, as it can be 
observed by considering the classical dependence of 
DNA melting temperature on pH and salt concentra- 
tion. In the present work, we report a single-molecule 
study of the stability of DNA against denaturation 
when subjected to changes in solvent. We investig- 
ated the appearance of DNA instability under specific 
external applied force and imposed twist values, 
which was revealed by an increase in the temporal 
fluctuations in the DNA extension. These fluctuations 
occur in the presence of a continuous interval of 
equilibrium states, ranging from a plectonemic 
state to a state characterized by denaturation bub- 
bles. In particular, we observe the fluctuations only 
around a characteristic force value. Moreover, this 
characteristic force is demonstrated to be notably 
sensitive to variations in the pH and ionic strength. 
Finally, an extension of a theoretical model of ple- 
ctoneme formation is used to estimate the average 
denaturation energy, which is found to be linearly 
correlated to the melting temperature of the DNA 
double strands. 



INTRODUCTION 

DNA denaturation is at the origin of many biological 
phenomena, for example, replication, transcription and 
the interaction of specific proteins with single-stranded 
DNA. In the present work, we conducted a direct quan- 
titative characterization of DNA denaturation at the 
single-molecule level. Given the influence of environmen- 
tal conditions on the biological activity of DNA, we per- 
formed measurements in different conditions of pH and 
ionic strength. 

Single-molecule techniques are particularly relevant in 
the current panorama of biophysical research because 
they have opened new roads to the analysis of several bio- 
logical and biomedical problems that are not easily observ- 
able with ensemble measurements (1-3). Typical examples 
of single molecule techniques are magnetic tweezers (MT) 
(4-8), optical tweezers (9-12), atomic force microscopy (13- 
16) and single-molecule fluorescence spectroscopy (17,18). 
Compared with the other techniques, MT methods have a 
straightforward implementation, and given their capability 
for applying torsional constrain, MT allows the study of 
single molecules in conditions that are not easily accessible 
with the other methods (19,20). Recently, valuable 
torsional experiments have also been performed with optical 
tweezers with unconventional nanoparticles (21-23) or by 
exploiting refined analysis (20,24), but MT nanomani- 
pulation still offers a simpler strategy for applying force 
and twist to single molecules. 

The MT technique is based on the manipulation of a 
single DNA molecule linked to magnetic spherical 
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microparticles by means of external magnetic fields (6). 
The microparticles are suspended in water solutions, but 
they are not freely diffusing because they are connected to 
the inner walls of a capillary tube by single DNA mol- 
ecules. Force and twist are applied to the bead through a 
field generated by external permanent magnets or electro- 
magnets. The movement of the bead is transferred to the 
sample, and thus force and twist are applied to the DNA. 

The first MT studies addressed the topology of individ- 
ual DNA molecules by measuring DNA extension as 
a function of the applied force and imposed torsion 
(25-27). In particular, it was observed that DNA exten- 
sion is influenced by torsion producing plectonemes (28). 
In certain conditions of torsional stress at high forces, MT 
nanonianipulation can also induce a denaturation, which 
is obtained only with the negative torsion because of the 
chirahty and natural hehcal orientation DNA (4). In the 
present work, we show the results of a series of MT ex- 
periments, where we evaluated the dependence of the sta- 
bihty of double-strand DNA, that is, denaturation 
bubbles (29-33), DNA breathing (34) or DNA melting 
(35,36), on the imposed turns and apphed force. The ex- 
periments were performed with DNA in standard phos- 
phate buffered sahne (PBS) at pH 7.4 and at various pH 
values and ionic strengths. The results are particularly 
interesting because DNA stability influences a host of bio- 
logical functions, such as DNA binding with transcription 
machinery and the replication process (37-39). It should 
be observed that the energy landscape of the denaturation- 
plectoneme transition was recently analysed in (40), and in 
the present article, we extend those results probing expli- 
citly the effects of pH and ionic strength. Furthermore, 
other articles discussed the DNA phase diagram in the 
force-torsion space (35). Here, we concentrate on the 
plectonemic-denaturation transition, in the force-turns 
space, never explored before. 

MATERIALS AND METHODS 

Capillary flow cell preparation 

The capillary flow cell (1x1 mm^ section, 5 cm long, 
VitroCom, Mountain Lakes, NJ, USA) was prepared ac- 
cording to a method in the literature (41). Briefly, the 
internal walls were uniformly coated with polystyrene 
(average molecular weight = 230000, Aldrich, Milan, 
Italy) (42) and were coated with 5 |.ig of sheep polyclonal 
anti-digoxigenin antibody (Roche, Milan, Italy). The 
treated surface was passivated for 2h at 37°C with a 
solution consisting of PBS pH 7.4, 0.1% of Tween-20, 
lOmg/ml of bovine serum albumin (Roche) and 3mM 
of NaNs (4). 

DNA and magnetic bead preparation 

The DNA molecule required to perform an MT experi- 
ment must be bound at one end to an immobile support 
(the glass capillary) and at the other end to a magnetic 
bead (5). The DNA molecule used in this work was 
7000 bp long and was assembled by a central part ex- 
tracted from pCMV6-Neo plasmid (5780 bp, ~1.9|im) 
and by two functionalized ends. One end was modified 



by digoxigenin, and the other end was modified with 
biotin to provide a connection to the immobile support 
(the anti-digoxigenin-coated capiUary) and to the mobile 
support (the avidin-coated magnetic bead), respectively. 
Following the methods previously reported in the litera- 
ture (6), the DNA tails were functionalized by polymerase 
chain reaction; thus, the DNA was attached via multiple 
bonds to the supports. The molecule results torsionally 
constrained so that the rotation of the magnets induces 
a real torque on the DNA, which could not swivel around 
the bonds. As it was already reported (6), we used l-|im 
diameter, streptavidin-coated, superparamagnetic beads 
(Dynabeads MyOne Streptavidin CI, Dynal, Invitrogen, 
Milan, Italy). We added 1 |il of the bead solution to 10^1 
of the functionalized DNA solution for a final concentra- 
tion of 1 X lO'' molecules of DNA per microhtres. After 
5 min of incubation, the DNA + microbeads solution was 
dfluted in 500 (il of PBS pH 7.4, supplemented with 0.1% 
of Tween-20 and 3 mM of NaN3 (referred to hereafter as 
PBS-Tween-20); the solution was subsequently injected 
into the capiUary. Finally, after 1 h of incubation, we 
washed out all of the unattached beads with a solution 
of PBS-Tween-20 pH 7.4. 

pH and ionic strength characterization 

To study the influence of environmental conditions, DNA 
properties were analysed in solutions having different pH 
and ionic strength values from those of the standard PBS 
at pH 7.4. The PBS solutions for measurement in different 
pH conditions were prepared by balancing Na2HP04 and 
KH2PO4 to obtain different values of pH and by changing 
the NaCl concentration to maintain constant ionic 
strength. The resulting values of pH were checked by a 
standard pHmeter (Mettler Toledo, Milan, Italy). The 
range of pH explored was 5 < pH < 9. 

The PBS solutions for measurement in different ionic 
strength conditions were prepared by modifying only the 
NaCl concentration (5, 50, 150, 395 and 500 mM) and 
maintaining the pH at 7.4. The resulting conductivity of 
the solution was checked by a conductimeter. For meas- 
urements at different conditions, at least 10 volumes of the 
solution was injected into the capillary to ensure complete 
substitution of the buffer. 



MT technique 

Our MT apparatus followed the most classical proposed 
schemes (5-7), and it has been previously described else- 
where (41,43,44). The capillary, which is hnked to a 
buffer-flow system, is placed over a homemade inverted 
microscope, which is composed of an objective (Nikon 
lOOx, 1.25-NA oil immersion with a 15-cm focal-length 
tube lens leading to an actual magnification of 75x) and is 
equipped with a piezoelectric focusing system (PIFoc, 
Physik Instruniente, Bresso, Italy). Images were acquired 
at a rate of 60 frames per second with a Charge-Coupled 
Device (CCD) camera (Marhn, Allied Vision Techno- 
logies) and were fed into a Personal Computer. 
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Data analysis 

By considering the diffraction images generated at differ- 
ent heights of the bead, illuminated by a hght emitting 
diode (LED) light (45,46), the DNA extension was 
calculated. A force is exerted (by the magnetic field of 
the external magnets) on the bead and consequently also 
on the DNA; the force was measured from the x-y fluc- 
tuations of the bead, as previously described (7,8). The 
mean square displacement a\ — {(x — {x))^) was evaluated 
with an appropriate correction for the camera integration 
time (8,41,47). Next, we extracted the force F with the 
equipartition theorem (4,7): F= kg ■ T ■ L^,/a\, where kg 
is the Boltzmann constant, T is the temperature and is 
the DNA extension. Once the applied force was calibrated 
in this manner, we were able to apply a known force to the 
bead. In the explored regime, the DNA extension L^. is an 
increasing function of the applied force F and the specific 
dependence of versus F can be interpreted in the frame- 
work of the so-called worm-like chain model (48). By 
simply rotating the external magnets at a fixed height 
(and thus at a fixed applied force), we can apply a twist 
to the torsionally constrained DNA molecule. 



RESULTS 

We have measured the average extension (L^.) of a DNA 
molecule as a function of the number of imposed turns n, 
or analogously as a function of the degree of supercoiling 
a — (LkQ is the natural linking number in absence of 
external constraints, and AL/, is the excess or deficit of 
L/,). The obtained data are shown in Figure 1 (upper 
panel) for the classical representative case of DNA at 
standard conditions (pH = 7.4, ionic strength 
Is = 1 50 mM) and for two different values of the applied 
pulhng force {F = 1.2 and 0.45 pN). The mean value (L^) 
was obtained by averaging over a time interval of several 
seconds, as detailed later. The circular labels A, B, C and 
D located on the data of Figure 1 (upper panel) corres- 
pond to the sketches of Figure lA-D (lower panels), 
which provide an idea of the actual behaviour of DNA 
under the actions of torsion and force. As previously 
described in literature (6) and visually depicted by the 
sketches of Figure 1, the effects of applied force and 
imposed twist are as follows: forming a plectoneme 
(Figure IC), extending the DNA (Figure lA and B) and 
forming a denaturation bubble (Figure ID). There are 
detailed experimental studies of the initial stages of 
plectoneme formation and rigorous nanomechanical 
models of the buckhng and post-buckhng behaviour of 
chains under torsion (49-53). However, we beheve that 
our naive sketches are still valuable and describe the 
primary characteristics of the experimental findings. 

The (Lp) as a function of «, data are presented in 
Figure 2 for certain values of the apphed pulhng force 
(F= 1.2 pN, F= 0.8 pN and F= 0.5 pN). The results 
correspond to those of the classical MT hterature (6). In 
this study, we concentrated our attention on the temporal 
statistical fluctuations of the data and particularly on the 
standard deviation ui^, of the DNA extension measure- 
ments. The histograms in the insets of Figure 2 show the 
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Figure 1. Upper panel: representative average DNA extension (L^,) 
measured using MT as a function of tlie number of imposed turns n, 
(lower X axis) and of the degree of supercoiling a (upper .t axis). Data 
obtained in standard conditions (pH = 7.4, = 150mM) with F= \ .2 
pN (upper data) and F = 0.45 pN (lower data). Lower panels A, B, C 
and D: schematic sketches of the DNA-magnetic bead arrangements 
during the MT experiments. The red circles over-plotted on the data 
highlight the correspondence between the schematic sketches and the 
DNA structures. 

distribution of ALg — (Lp — {L^)) for the data highlighted 
by the coloured circles. The histograms are fitted by 
Gaussian distributions, and the resulting values of the 
standard deviations a^^. for F = 0.8 pN and n, < —20 
(i.e. in the transition region between the plectonemic and 
denaturation states) seem to be larger with respect to other 
domains of F and 

To explore this phenomenon quantitatively, we also 
calculated the temporal autocorrelation function 
g(r) — {Le(t)Le(t+r)) / {Li.(t)') for various values of the 
time delay i. Two examples of the obtained correlation 
functions and the corresponding original data traces 
of the DNA extension versus time are presented in 
Figure 3, panel A and B, respectively. As apparent from 
Figure 3, the correlation function is negligible for F = 0.5 
pN, that is, the correlation between the data occurs on a 
short time interval compared with our acquisition time. 
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Figure 2. Average DNA extension (L^) measured using the MT technique as a function of the nuinber of imposed turns n, (lower x axis) and of the 
degree of supercoiling a (upper .t axis) at the various forces F indicated by the labels. The lateral histograms represent the statistical distribution of 
ALj = — (L;,) (in |im) of the data highlighted by the coloured circles. Data obtained in standard conditions. 
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Figure 3. Panel A: normalized autocorrelation functions 
g(T) = {Lg(t)L,,(t+T))/{L,,(t)-) of the DNA extension fluctuations 
measured as a function of the delay time i (symbols) and exponential 
fitting of the data (continuous line), which allows extracting the char- 
acteristic Tj./,„r. Panel B: original temporal traces of the DNA extension 
fluctuations AL^, corresponding to the data of Panel A. Data obtained 
in standard conditions at «, = —40 and F ^ 0.5 pN (black symbols) and 
0.8 pN (red symbols). 



For F = 0.8 pN, the correlation function is a decreasing 
function over an appreciable duration time of the order of 
100 ms. We fitted the correlation functions with an 
exponentially decaying expression of the form 
^(r) = e~^, where !,./,„,. (the characteristic correlation 
time) is the free fitting parameter. 

The data of the standard deviations obtained as a 
function of the number of imposed turns and the 



applied forces F are presented in Figures 4A and 5A, 
respectively. Analogously, the measured !,./,„,. values, as 
resulting from the previously described fitting procedure, 
are shown in Figures 4B and 5B for different values of the 
Ht and F. In Figure 4A and B, we have also showed the 
error bars as coloured bands around the data. Given 
the acquisition frame rate of the camera, we choose a 
value of 50 ms as a confidence threshold. The data below 
that threshold are not reliable. To emphasize this conside- 
ration, we have added to the figures showing !,./,„,. a dashed 
horizontal fine, indicating the limit of validity of the 
experimental fitting procedure for extracting t^/,„,.. 

It is worthwhile to note that to avoid systematic statis- 
tical errors, the values of cr^^, and {Lg) are obtained by 
making a statistical average, that is, collecting data 
during a time interval of ~100 of the corresponding 
T,./,„,.. However, we have verified that the plotted values 
of CTi^, T^./,„,. and (Lp) were not modified by increasing 
further the time interval of the statistical averaging. As 
clearly shown by Figures 4 and 5, the data seem to be 
more noisy and more correlated in a specific range of 
and F. Indeed, the values of the standard deviation ctl^ 
and T,./,„,. below a fixed negative value of n, (rit < —20) 
are larger with respect to other domains of F and 
The values of aw and !,./,„,. show both a clear maximum 
around a characteristic value of the force = 0.8 pN). 
The value of i^c/iar coincides for both a^^, and Tj./,„,.. The 
obtained F,./,,,^ seems to not depend on «, within the ex- 
perimental errors, and the maximum values of a^p and 
Tp/,„,. at = F,./,,,,., in the range explored, are decreasing 
functions of 

The domains with large values of ct/.^, and Tj./,„,. are 
located in regions where, under the action of the 
imposed turns and forces, the DNA is fluctuating 
between the denaturated and plectonemic states, continu- 
ously exploring all the intermediate states. Because the 
possibility for DNA denaturation is strongly dependent 
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Figure 4. Standard deviation a^^ of the DNA extension (Panel A) and 
cliaracteristic time t^./,,,,. of the correlation functions of the DNA exten- 
sion fluctuations (Panel B) measured as a function of the number of 
imposed turns n, (lower x axis) and of the degree of supercoiling a 
(upper X axis). The horizontal dashed line represents the reliability 
threshold of the fitting procedure for obtaining T,./,„r, because of the 
finite acquisition time of the experiments. The coloured bands represent 
the error bars (confidence regions) as obtained by considering the 
quality of the fitting procedure. Data obtained in standard conditions 
at the different forces F indicated by the labels. 



on the acid-base or ionic characteristic of the solvent 
in which the DNA is dispersed (54), we have repeated 
the aforeinentioned analysis for different values of pH 
and I,. 

To investigate this point, we measured (L^) versus «, at 
fixed force for various values of pH and I,, and the results 
are presented in Figure 6A and B, respectively. In this 
study, we directed our attention to the data at negative 
rit, where a mild or flat dependence of {L,,) as a function of 
\n,\ is a clear indication of the induced mechanical de- 
naturation of the DNA molecule. As noticeably 
demonstrated in Figure 6A, the denaturation occurred 
more easily for pH values far from pH 7.4, where the 
hydrogen bonds of the bases are weaker (55). Similarly 
shown in Figure 6B, the DNA is easily denaturated for 
negative at 1, = 5 mM, that is, at low ionic strength, 
where the screening of the electrostatic repulsion of the 
two double strands is reduced. 
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Figure 5. Standard deviation ai,^ of the DNA extension (Panel A) and 
characteristic time Xci,ar of the correlation functions of the DNA exten- 
sion fluctuations (Panel B) measured as a function of the force F. In all 
cases, the estimated errors are below the 20% of the measured values. 
The horizontal dashed line in Panel B represents the reliability thresh- 
old of the fitting procedure for obtaining T^itan because of the finite 
acquisition time of the experiments. Data obtained in standard condi- 
tions at the different number of imposed turns n, indicated by the 
labels. The vertical line indicates the characteristic force -F,./,„r. 



More quantitatively, in Figure 7A-D, we have plotted 
the extracted values of the standard deviations and the 
characteristic correlation times measured as a function 
of the applied F for different pH and 1, at negative 
imposed turn n, = —40. From the entire ensemble of 
data, we extracted the values of the characteristic force 
in the different situations. The data show that at pH 
values other than pH 7.4 or for low ionic strength, the 
DNA transition to denaturation is shifted to lower 
values of i^(./,„,-, that is, the DNA is less stable. By consider- 
ing Figure 7, we can observe that the distribution of forces 
over which the DNA length fluctuations occur is system- 
atically larger for high ionic strength (and less systematic- 
ally for non-physiological pH). Indeed, the denaturation 
occurs at lower forces in acid or basic environments or at 
low ionic strength, as shown in Figure 8A and B, where 
the F,./,ar values are plotted as functions of pH and 1,. 

The observed bell-like shape of F,.,,^,- versus pH and the 
increasing dependence of F^./,,,,. versus I, are notably similar 
to the classical dependences observed for the DNA 
melting temperature T,„ as functions of pH and 
Indeed, in Figure 8C and D, for comparison we plotted 
the behaviour of T,„ versus pH and 1, as classically 
reported in hterature (56,57). The similarity between the 
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Figure 6. Average DNA extension (L,,) measured as a function of the 
number of imposed turns ii, (lower x axis) and of the degree of super- 
coiling a (upper x axis). Data obtained at different pHs (see labels) 
with a constant force f =0.75 pN and ionic strength /, = ISOmM 
(Panel A). Data obtained at different ionic strengths /, (see labels) 
with a constant force F= 0.8 pN and pH = 7.4 (Panel B). 



behaviours confirms the capabihty of the method in 
characterizing the stability of DNA in different situations. 



DISCUSSION 

In the experimental results presented in this report, we 
mainly described the behaviour of DNA for negative 
imposed turns and for characteristic values of the 
appHed force, where the DNA extension fluctuations are 
largely enhanced; the enhanced fluctuations are observed 
because molecules can choose between the mechanical 
options of locally denaturing (bubbles) or forming super- 
coils (plectonemes). We observed that such fluctuations 
are strongly dependent on the ionic strength and on the 
pH of the solvent. 

We now recall the main results of a simplified quanti- 
tative model introduced in (40) which can describe all of 
the experimental findings. According to that model (40), 
when the DNA is twisted, the phenomenon of denatur- 
ation or plectoneme formation is accompanied by the 
stored energy .E,,,, = £:„„,r+£'p/m+£rf™ = 5£(27r ■ «<,)^+ 

{InRp ■ tip) ■ (2I2 +F)+a ■ n^. Where the first term is the 



twisting energy (£■,„.,-.„), the second is the energy 
accumulated during plectoneme formation {Epi^^,) and 
the third one is the denaturating energy Lq is the 

extended length of the DNA, — {tit — np — rid) is 
the number of turns that store energy in the twist, rip is 
the number of loops of radius Rp along the DNA, nj is the 
number of turns relaxed by a partial DNA denaturation, 
B is the DNA bending constant, C is the DNA torsion 
constant and the constant a is proportional to the binding 
energy between the bases of the two DNA single strands. 
The first two terms were already used in model (6), 
whereas the last one was originally introduced in (40). 

The plectonemic energy has a double origin: the work 
necessary for bending a lit ■ Rp portion of the double 
strand and the work for puUing the bead in presence of 
the force F. Disregarding the denaturation contribution, 
by minimizing the energy E = E„rjst+Epif,ct, one can obtain 
the standard results of the equihbrium plectoneme radius 

Rp — and the buckling number «/, = Lq ^^c^' , as pre- 
viously reported in the literature (6). 

The situation is more comphcated when the denatur- 
ation contribution energy E,ie„ is explicitly taken into 
account. The boundary conditions prescribe for E,o, a tri- 
angular validity region in the plane {np,n,/) with 
rit < rip < 0, nt < rid < 0 and rip+Hd < The model (40) 
predicts the presence of a characteristic force E,./,^,., where 
the energy necessary for relaxing a new turn by forming a 
loop or by denaturing the double strands is equal. The 
model analytically relates this characteristic force to the 
bending constant B and to the parameter a: Fcimr = gf^- 

The representative energy landscapes for three different 
values of the external force E are plotted in Figure 9A 
(E < Echar), Figure 9B (E = E,i,„r) and Figure 9C 
(E > Echar) as a function of iip and «,/. The darker 
regions in Figure 9A, B and C represent the equihbrium 
points or hues of minimum energy. The energy landscapes 
are obtained assuming standard values of 
5= 50■^'B■ r,C= 140■^'B■ r, Lo=1.88 ^m and 
a — 1.12- 10^'^J. In the first case, when E < Echar. the 
energy landscape presents a minimum for iip — n, — 1% 
(i.e. plectoneme formation) and «,/ = 0 (i.e. no denatur- 
ation). In the case when E > E^har, the minimum of the 
energy is attained for — n, — Lq ■ a/{ACjp-) (i.e. the 
force is enough to form denaturation bubbles) and 
Up = 0 (i.e. no formation of plectonemes). In the inter- 
mediate situation, when E — Eduir, the minimum of the 
energy landscape is represented by a line, suggesting that 
the force is allowing a fluctuation along the line connect- 
ing the two states. 

Using the experimental value E^ha,- ~ 0.8 pN and a 
DNA persistence length of ~50nm, the resulting value 
of a is ~1.1 ■ 10~'^ J 66kJ/mole), which is appro- 
ximately eight times the classically reported average 
value of the free energy per DNA base pair, AG ^ 8.4 
kJ/mole (58-60). It seems then necessary to open eight 
base pairs to relax one turn. This number of base pairs 
is consistent with the pitch of the double helix. This value 
is obtained not considering the influence of other possible 
structures [e.g. hairpin (61) or cruciform (62)]. The 
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Figure 7. Standard deviation at, of the DNA extension (Panels A and B) and cliaracteristic time Tci,ar of tlie correlation functions of the DNA extension 
fluctuations (Panels C and D) measured at the different values of pH (/, = 150mM) and ionic strengths (pH = 7.4) indicated by the labels. All data are 
acquired at negative imposed turn n, = —40, and in all cases, the estimated errors are below the 20% of the measured values. The horizontal dashed line 
represents the reliability threshold of the fitting procedure for obtaining x^./ja,-, because of the finite acquisition time of the experiments. 
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Figure 8. Characteristic force frf,„r measured as a function of the pH 
(Panel A) and of the ionic strength (Panel B). Literature data (56,57) 
of the DNA melting temperature T„, measured as a function of the pH 
(Panel C) and of the ionic strength /, (Panel D). 



obtained value of a can be easily compared with the 
measured denaturation torque reported in (63). The par- 
ameter a is obtained by calculating the work necessary for 
twisting the DNA of one turn: a —T ■ In (where F is the 



torque). By considering the reported torque of llpN- nm, 
the calculated value for a is 0.7-10~'V, consistent with 
what observed in the present work. 

Using the proposed energy expression E,o, and a 
Boltzmann distribution enables a prediction for {L^) and 
a^^, in the different situations. The calculated average 
DNA extension {Lf,) is presented in Figure 9D as a 
function of n, for three different force values. In this 
study, we suppose that the DNA extension variations 
are mainly because of plectoneme formation 
(Lp = Lo — IttR ■ rip). We also assume that the denatur- 
ation phenomenon does not introduced a significant 
DNA length variation given the low n, values at which 
the experiments are performed. These n, values correspond 
to a ssDNA of 300-350 bp, which is <5% of the DNA 
length. Evaluating the extension length of such ssDNA, 
and assuming 1 nm of persistence length, the maximum 
extension variation should be of the order of few 10th of 
nanometres, which is not relevant in our case. In the 
present calculations, we also introduced a worm-like 
chain model dependence of the extended value of the 
DNA length Lq (64), and we also verified that the force 
dependence of the constant C is not strongly influencing 
the results (65). Comparing Figures 2 and 9D, we note 
that the model qualitatively describes not only the classical 
plectonemic behaviour observed for n, >0 but also the 
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Figure 9. Contour plot of the theoretical energy landscape simulated as a function of the number of turns np relaxed in plectonemes and the number 
of denaturated turns iij, calculated with the expression for E,„, described in the text, assuming 0.7 pN = F < F^./,,,, (Panel A), 0.8 pN = F = F^.i,,,, 
(Panel B) and 0.9 pN = F > F^.i,,,,- (Panel C). Intensity colour scale is expressed in terms of kg ■ T. Simulations obtained with the following param- 
eters: B = 50 ■ ks ■ T,C = 140 ■ ks ■T,Lo = 1.88/im, a? = 1.12 ■ 10""J. Panel D: calculated average DNA extension (i,,); Panels E and F: calculated 
standard deviations a^^ of the DNA fluctuations. Data calculated with the model described in the text as a function of the number of imposed turns 
n, for the different values of the applied force F indicated in the labels. 



transition between the plectonemic and denaturation 
regions for n, <0. 

The quahtative agreement between the data and the 
model is also demonstrated by comparing the behaviour 
of the calculated cr/. , as a function of «, shown in Figure 9E 
and the experimental findings shown in Figure 4A. Both 
the calculated and measured fluctuations increase for 
«, <0 and in a specific range of F. The qualitative agree- 
ment is confirmed by Figures 9F and 5A, where we have 
shown the calculated and measured values of ct^,, as a 
function of F. The similarity of the theoretical predictions 
and the data in all of the situations confirm the capability 
of the model to describe the characteristics of the experi- 
mental findings. 

There are certain obvious quantitative differences 
between the measurements and the theory. For example, 
the calculated fluctuations show a transition from the 
plectonemic state towards the denaturating state in a 
narrow transition region of the force (±10% of 
the i^i /,„,), whereas the experiments report a range larger 
(± 20% of the F,./,^,) than the model. This larger experi- 
mental fluctuation amplitudes with respect to the theoret- 
ical ones are presumably because of the fact that our model 
neglects several aspects of the real system, such as the 
presence of an initiation energy for starting the denatur- 
ation bubble (52), the possibility to have different bubbles 
along the DNA, the occurrence of different plectonemic 



regions (66) and the existence of thermal oscillations of 
the plectoneme radius (53). These considerations are also 
presumably the reason why the proposed model does not 
explain the observed dependence of the distribution of 
forces, over which the DNA fluctuations occur, on ionic 
strength or pH (Figure 7). In particular, an influence of the 
buffer condition on the number of plectoneme domains has 
been recently demonstrated (67). All these not considered 
ingredients can be at the origin of the quantitative differ- 
ences between data and simulations. Nevertheless, the pre- 
sented results reproduce qualitatively all the experimental 
conclusions, and the modelling can provide a deeper insight 
into the characteristics of DNA denaturation under the 
effects of external constraints. Other DNA structures 
(Z-DNA or L-DNA) can be formed at n, < 0 (68-70), 
but all such structures are present in non-physiological con- 
ditions, far from the low force and turn regime, which is the 
object of the present work. 

To consider the data taken in non-standard conditions of 
the solvent, we recall that the classical phenomenon of 
DNA denaturation that occurs at the temperature T,„ 
depends on the specific DNA base composition and on 
environmental conditions (54). Indeed, at the temperature 
T,„, the hydrogen bonding between the bases weakens, and 
the two strands separate. T,„ is a strongly dependent func- 
tion of the pH and the ionic strength, as reported in the 
classical measurements (57,56) shown in Figure 8C and D. 
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Figure 10. Symbols: measured denaturating parameter ot plotted as a 
function of the melting temperature T„, obtained at different pH and 
ionic strengths /, from the literature (56,57). Continuous line: linear 
fitting of the data. The error bars correspond to the distribution 
width of the a values obtained probing different molecules (~10). 



We assume that the DNA mehing occurs when the 
number of bases in the denaturation bubbles is compar- 
able with half of the number of DNA base pairs. 
Considering a simplified version of the Poland-Scheraga 
model as reported in (71), the melting temperature T„, is 
related to the energy E^,/ for elongating a bubble by one 
base pair, by the relation = 2 ■ ln{2). Considering 
Eel — oi/%, where 8 is the numl^'er of base pairs opened 
to relax one loop, we predict a linear relationship 
between a and T,,, with a theoretical proportionality 
constant m of ~1.5 ■ \0~^^J/°K. In Figure 10, the values 
of a, which are plotted as a function of the melting tem- 
perature measured at different pH and confirm this hn- 
ear behaviour with a constant m = 3.4± 1 ■ \Q^^^J/°K. 
The order of magnitude agreement validates the appropri- 
ateness of the proposed model. 



CONCLUSION 

The data and the corresponding model shown in the 
present work elucidate new, innovative and fundamental 
aspects of DNA biophysics, specifically, the nano- 
niechanical stabihty of single DNA molecules under 
various environmental conditions in a force and torsion 
regime compatible with those involved in important 
protein activities, such as DNA and RNA polymerases 
(72,73) and DNA gyrase (74). The presented research 
opens the way for further studies and investigations in 
more sophisticated and realistic situations; these studies 
could examine modified DNA structures or the effects of 
the presence of DNA binding proteins, enzymes and 
Ugands. 
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